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Abstract 

We investigate a curvature singularity appearing in the star collapse process in F[R) gravity. In 
particular, we propose an understanding of the mechanism to produce the curvature singularity. 
Moreover, we explicitly demonstrate that ii" (1 < a < 2) term addition could cure the curvature 
singularity and viable F{R) gravity models could become free of such a singularity. Furthermore, 
we discuss the realization process of the curvature singularity and estimate the time scale of its 
appearance. For exponential gravity, it is shown that in case of the star collapse, the time scale is 
much shorter than the age of the universe, whereas in cosmological circumstances, it is as long as 
the cosmological time. 
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I. INTRODUCTION 



Modified gravity is expected to be the fundamental scenario to give tlie unified gravita- 
tional alternative to dark energy and inflation A number of viable modified gravities 
are known. Of course, there exist important criteria for viability, such as the fulfillment of 
Solar System tests. Among these criteria, one of the most important ones is related with 
so-called matter instability [2] in F{R) gravity, which means that the curvature inside mat- 
ter sphere becomes very large, i.e., strong gravity (for an introduction to F{R) gravity and 
very recent reviews on it, see, e.g., [l, 3|). It was indicated that such matter instability may 
be dangerous in the relativistic star formation process {4] due to the appearance of corre- 
sponding singularity. On the other hand, the hydrostatic equilibrium of a stellar structure 
in the framework of F{R) gravity has been explored by studying the modified Lane-Emden 
equation in Ref. jsl. 

Recently, the instability in F{R) gravity has been discussed for a gravitating system with 
a time dependent mass density like astronomical massive objects in Ref. js]. It has been 
shown that when a star shrinks and the mass density becomes larger, the time-dependent 
matter instability, in which the scalar curvature becomes very large, could occur for a class 
of F{R) gravity models. It is interesting to understand how common such time-dependent 
matter instability is in F{R) gravity and how viable F{R) gravity models may be protected 
against it. In this paper, we study the generation mechanism of the time-dependent matter 
instability in the star collapse. We show that the time- dependent matter instability develops 
and consequently the curvature singularity could appear in a viable F{R) gravity model [7| 
and some version of exponential gravity SMlOj . We note that the equivalent or very similar 
modification of gravity to Ref. Q has been considered in Refs. [l2|. In addition, we 
demonstrate that the curvature singularity could be cured by adding the higher derivative 
term i?" (1 < a < 2) and viable F{R) gravity models could become free of such a singularity. 
Furthermore, we discuss the realization process of the curvature singularity in the viable 
F{R) gravity model and estimate the time scale of its appearance in exponential gravity. 
We show that in case of the star collapse, the time scale is much shorter than the age of 
the universe, whereas in cosmological circumstances, it is as long as the cosmological time. 
We mention that the problem of singularity in the star collapse was studied in Ref. jl^, in 
which it was stated that the curvature singularity in the future does not appear if the model 
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of F{R) gravity is built very carefully (for an earlier proposal, see [ij]). It is considered that 
the curvature singularity could emerge in a generic F{R) gravity model unless fine-tuning 
is taken. We use units of k-Q = c = h = 1 and denote the gravitational constant SnG by 

= Svr/Mpi^ with the Planck mass of Mpi = G'^/^ = 1.2 x lO^^GeV. 

The paper is organized as follows. In Sec. II, we first review the matter instability and 
study the generation mechanism of the time-dependent matter instability. In Sec. Ill, we 
examine how the curvature singularity occurs and analyze the time scale of its appearance. 
Finally, conclusions are given in Sec. IV. 



II. CURVATURE SINGULARITY IN THE COLLAPSE OF STARS 
A. Matter instability 



We start with reviewing the matter instability issue [2| in F(R) gravity. It is related 
with the fact that spherical body solution in general relativity may not be the solution in 
modified theory. It may appear when the energy density or the curvature is large compared 
with the average one in the universe, as is the case inside of a star. 

We consider the following action 

S = j rf^Xv/^ 1^ + £^attcr| , (2.1) 

where g is the determinant of the metric tensor gf^j,, F{R) is an arbitrary function of R, and 
£matter IS the matter Lagrangian. The trace of the gravitational field equation derived from 
the action in Eq. fl2 . 1 p is given by 



F(3)(i?) F'{R)R 2F{R) 



K 
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+ 3F(2) (R) ^P^^'^ + 3^7^(2) (^) 3^(2) (^) - 3^(2) (^) ^-atter , (2.2) 

where a prime denotes a derivative with respect to R, is the covariant derivative operator 
associated with g^,^, □ = g'^'^'V^Vu is the covariant d'Alembertian for a scalar field, and 
Tjnattcr = ^matter p'' IS the tracc of the matter energy-momentum tensor T^a^tter fiu defined as 
Tn^atter;.. = ' (2 / (5 C„,,ttcr / 5 g^ ■ We also denote d^F{R)/dR^ by F^^\R). Let us 
now examine the perturbation from the solution of general relativity. We express the scalar 
curvature solution given by the matter density in the Einstein gravity by R\, ~ K^pmatter > 
with pmatter being the energy density of matter and separate the scalar curvature R into 
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the sum of i?b and the perturbed part -Rp as i? = -Rb + -Rp (I-Rpi - C [-R hD- Substituting 



this expression into Eq. (12.21) . we obtain the perturbed equation 15|, |l6|. It is convenient 
to consider the case that -Rb and -Rp are homogeneous and isotropic, that is, they do not 
depend on the spatial coordinates. Hence, the d'Alembertian can be replaced with the 
second derivative with respect to the time coordinate: D-Rp — —d^Rp. As a result, the 
perturbed equation has the following structure: 

= -d^Rp + U{Rh)Rp + const. . (2.3) 

Thus, if W(-Rb) > 0, -Rp becomes exponentially large with time t'. Rp ~ 
system is unstable. 

For example, in the model of F{R) = R — a/R + bR^ with a and b being positive 



constants 



15| . for b ^ o/l-Rbl one gets 



W(i?b)~^>0. (2.4) 



Therefore, the system could be unstable. Since -Rb is estimated as [2|, however, 

R^ (lOhec)-' ( , (2.5) 
^ ^ \gcm~'^/ 



a 



the decay time is ~ 1, 000 sec, i.e., macroscopic. On the other hand, in a viable model 

FiR)=R-m^ (2-6) 
C2 [R/m^) + 1 

where Ci and C2 are dimensionless parameters, n{> 0) is a positive constant, and m is a mass 
scale, if one takes ci > 0, U{Rh) is negative as 

m^.)-- , <0- (2.7) 

n[n + l)ci 

Consequently, the system could be stable and there is no matter instability. 

We remark that — W(-Rb) can be regarded as the square of the effective mass uiIq for 
the scalar mode -Rp, given by m^g ~ F"{R)/3 [l7]. This means that if F"{R) > 0, one 
has m^fj > and hence the scalar mode -Rp is stable, whereas it is unstable for F"{R) < 
because m^g < 0. If F"{R) < 0, this instability occurs inside matter when the deviation of 
-R from -Rb ~ /t^Pmatter emerges. Thus, it can be interpreted as a matter instability. 
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B. Time-dependent matter instability and curvature singularity 

In Ref. lal, it has been shown that even for the viable models such as those in Refs. 
the increasing matter density could generate the curvature singularity, which we will describe 
in this subsection. We write F{R) as^ 

F{R) = R + f{R). (2.8) 

The trace equation of the gravitational field is described as 

SDf (i?) -R + Rf'{R) - 2f{R) = ^'T^^tter • (2.9) 

We now consider a small region inside the star, where the energy density can be regarded as 
homogeneous and isotropic as assumed in Ref. I2I. Since we study the gravitational field of 
the star such as the Sun, whose mass density is p© = 1.4gcm~'^, we also suppose that the 
gravitational field is weak and the curvature is small because Rh ~ ^^Pmatter- (If a neutron 
star has the solar mass Mq = 2.0 X 10^3 g and the radius 1.0 x 10*^ cm, the mass density is 
4.7 X lO^^gcm"^, which is much larger than the solar mass density. The gravitational field 
of a neutron star is regarded as strong because a general relativistic treatment is necessary 
to explore a neutron star.) Then, we may replace □ in Eq. (12.91) with —df. By defining ip 
as (f = —f'{R), we can (at least locally) solve R with respect to ip: R = R{ip). Using this 
solution, Eq. (12.91) can be rewritten as 

^ = J^{ip,t) = ^{Riip) + R{ip)ifi + 2fiRiip)) + K^T^,ttcrit)} , (2.10) 

where the dot denotes the time derivative of d/dt. Compared with Newton's equation of 
motion mx = F, we can regard J^{(p,t) as a "force". From Eq. (I2.10p . we also obtain 

O = i02_ /■ dtmnm^t). (2.11) 



We explore the case in which Tmatter increases with time: 

t 

to 



^matter — ~2^matterO I ^ + ^j" ) ! (2-12) 



where Tmattero is constant and tn is a time. 



^ Note that the definition of F{R) and f{R) are different from those in Ref. 
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.BQ, 



In case of the viable models, e.g., those in Refs. [7|, f{R) is given by 

f{R)^-f, + A^, (2.13) 
nit" 

for large R. Here, /o and /i are constants. If f{R) = F{R) — R is written by Eq. (12. 6p . one 
finds /o = (ci/c2)m^ and fi/n = (ci / C2) m'^^^'^^K In this case, we have 

^ = ^' (2-14) 
and therefore J^{f, t) in Eq. (I2.10p is given by 



Combining Eqs. (12.1 ip and (I2.15p . we obtain 



(2.15) 



+ 



fi^^^^matterO ( 1 + ^ ) + 2/o 



^_ ^^T^atterO ^ ^^^^^^ ^ ^2.16) 



where £^ is a constant corresponding to "energy" in the classical mechanics and U{(p) to 
the potential. The case that R is large corresponds to the case that (p is small. Hence, we 
neglect the second term on the right-hand side of the expression for U {(p,t) in (I2.16p . i.e., 
- (1/3) [ff{n + 1)/ {2n + 1)] (1 + 2/n) ((^//i)^^"+^)/^"+^^ First, we consider the case that 
^matter = 0. It follows from Eq. (I2.16P that the "potential" U{!f) is given by 

C(^) = l/_!i±if^)=^+2/„4, (2.17) 



3 1 n \fi 
which vanishes at 

and has a minimum at 

if = ^ /r" (2/o)-("+^) . (2.19) 

The conceptual form of U (y?) is shown in Fig. [T] It is clear from Fig. [T]that if we start with 
an initial condition, for example, (p > ipo and <f < 0, ip reaches = 0, which corresponds to 
the curvature singularity R = oo. 








FIG. 1: The "potential" U (ip) in Eq. ()2.17p . The vertical and horizontal axes show U (ip) and tp, 
respectively. 



FIG. 2: The change of the "potential" U (p) in Eq. (I2.16P with time t. Legend is the same as 



When Tmattero 7^ 0, the potential shifts up with time as in Fig. [2l By growing the energy 
density, the singularity R = oo can be realized easily. Even if Tmatter = 0, there could appear 
the curvature singularity but the increase of Tmatter, which may occur by the star collapse, 
becomes a trigger to generate the singularity. 

A prescription to avoid the singularity is to grow the potential at i? = oo. In order to 
realize the growing up of the potential, we add the following term, say, to f{R) in Eq. (I2.13p 








Fig.m 



f{R) ^ f{R) - hR 
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Here, we assume /2 is a positive constant: /2 > and a is a constant. If we choose a > 1, 
the added term dominates when the curvature R is large and we find 

if ~ af2R"-^ . (2.21) 

Hence, the case R ^ oo corresponds to the case ip — )■ +oo. In this case, J^{(p, t) in Eq. f l2.10p 
is given by 

a 

J-(^,t)~/2(a-2)f^V" . (2.22) 



Since a > 0, the magnitude of the "force" J-'{ip,t) becomes infinitely large when R — )■ +oo. 
When a < 2, the force is negative and works to decrease ip or the scalar curvature, and 
therefore the curvature singularity is not realized. In fact, from Eq. (12.111) we have the 
following form: 



1 -2 , „, > fia(a-l)(2-a) { <p 



f = -^' + f/(^.«). UW)^ " '^^j; ^^^j . (2.23) 

Thus, the "potential" U{ip) is positive as long as 1 < a < 2 and becomes infinite when 
V9 — )■ oo, which means that we need infinite "work" to arrive at oo. On the other hand, 
when a > 2, the force is positive and therefore the force makes the curvature infinite and 
the curvature singularity is easily realized. Note that even if we choose /2 to be negative, 
the result is not changed because the curvature singularity corresponds to (/? — ?■ — oo in this 
case although the sign of the "force" J^{ip, t) is changed. 

It is also important to remark the case of adding R^ term, i.e., a = 2 as f{R) — )■ 
f[R) + /s-R^, where /s is assumed to be a positive constant. For the "realistic" model in 
Eq. (I2.13p . since the inverse power term /i/ {nR^) can be negligible when the curvature 
is large, we may consider the model F{R) R + f^R^- In this case, we may regard 
if ~ —2/3/2. In the expression of the "force" J-'(v5, t) in Eq. f l2.10p . the second term R {ip) ip 
is cancelled by the third term 2/ {R (v^)), which is consistent with Eq. (I2.22p . and the fourth 
term K^Tmatter(^) could be neglected. Hence, only the first term R^^p) contributes to the 
dynamics and it gives the "force" linear to (p as J^{(p,t) ~ —ip/ (6/3). Thus, from Eq. (12.1 ip 
we obtain quadratic "potential" U{(p) ~ ip'^/ (I2/3), which becomes positive infinity when R 
goes to positive infinity {ip goes to negative infinity). As a result, the curvature singularity 
is prevented. 
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Instead of the model in Eq. fl2.13p . we may investigate the exponential gravity model 

Q, 

/(i?) = -/,(l-e-il) , (2.24) 
where fc and Rc are positive constants. In the model in Eq. f l2.24p . we have 

if = ^e~^ . (2.25) 

Therefore, the curvature singularity R — )• +00 corresponds to (f = +0. For the model in 
Eq. (1221, we find 

^(^,t)~-f In^. (2.26) 

When {p — +0, the "force" J^{(p, t) becomes positive and infinite, but now it follows from 
Eq. (12. lip that we have the following form: 

S = \^' + UM, t/(^)~f ^(in^-l) , (2.27) 

which gives the finite value of "potential" when (f — )■ +0. This implies that we only need 
finite "work" to arrive at (f = +0. Therefore, the curvature singularity i? — i- 00 can be 
realized. This situation is almost the same as that in the model in Eq. f l2.13p because when 
(f — )■ +0, the "force" in Eq. f l2.15p diverges but the "potential" in Eq. f l2.16p is finite. 

We mention that there exists the following model with two exponential terms to realize 
inflation as well as the late-time cosmic acceleration 
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/(i?) = -/,(l-e-^) -/i 



1 - e (fli. 



+ -fRf^ , (2.28) 



where /i, -Ri(^ -Rc) and 7 are positive constants, /3 is a constant, and g(> 1) is a natural 
number. Here, we take -Rc and -Ri as the current curvature and the value of -R during 
inflation. If (i) R < R^ (ii) -Rc ^ -R ^ -Ri, and (iii) R ^ Ri, the first, second, and third 
terms on the right-hand side of Eq. f l2.28p becomes dominant, respectively. For the case 
(i), the model in Eq. f l2.28p corresponds to the one in Eq. f l2.24p . whereas for the case (iii) 
the behavior of the model in Eq. (12.280 is similar to that in the case of Eq. (I2.20p with 
the additional term domination. In terms of the case (ii), we examine the behavior for the 
corresponding form of f{R) as f{R) = — /i [l — e~^^^^^'^''~^ . In this case, we obtain 



where we have used R <^ Ri. For g > 1, increases as R becomes large. Moreover, we find 

The sign of the "force" J^(v5, t) is negative and therefore the force works to decrease 93. 
Furthermore, Eq. (12. lip becomes 

Since g > 1, f/ (</)) increases as (p becomes large. Consequently, the behavior of the model 
in Eq. (I2.28P in the intermediate regime /2c -R -Ri is stable. Note that this does not 
mean that the curvature singularity could not be generated because we are assuming the 
scalar curvature is finite as R <^ Ri although the curvature has a tendency to decrease by 
the "force" J'{ip,t). Without the last term in Eq. (I2.28p . there could appear the curvature 
singularity when R^ Ri. 

In Ref. [g], the emergence of the curvature singularity in the future has been investigated 
by the analytical studies, which have also been supplemented with numerical solutions. 
The equation for the evolution of R was reduced to the form of an oscillator equation 
X + dV{x)/dx = 0, where x and V corresponds to (f and U in this paper. In case of a time 
dependent potential V = V{x,t), since it is impossible to obtain an analytical solution, the 
qualitative behavior of the solution has been analyzed. 

III. REALIZATION OF THE CURVATURE SINGULARITY 
A. Realization process of the curvature singularity 

We discuss how the curvature singularity is realized. As an example, we investigate the 
viable model in Eq. (I2.13p . We consider a small region inside the star, which can be regarded 
as homogeneous and isotropic. In this case, the space-time is locally described by the flat 
Friedmann-Lemaitre- Robertson- Walker (FLRW) : 

ds^ = -df + a\t) J2 (^^')' ' (3-1) 

1=1,2,3 

where a{t) is the scale factor. Note that since we are examining the star collapse, the Hubble 
rate H = a/a is negative, that is, the space-time is shrinking. We also mention that the 
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energy densities of the matters automatically increase because the region is shrinking. Thus, 
we do not give an explicit time dependence as in Eq. fl2.12p . 

In case of cosmology, it is well-known that for the model in Eq. f l2.13p . Type II singularity 
can be realized, where the scale factor a and the effective energy density are finite but the 
effective pressure and the scalar curvature R diverge 18|. We apply this fact to the star 
collapse. Since a is finite, the energy density and the pressure from the matter are finite 
and can be neglected near the singularity. In such a situation, we find the Hubble rate H is 
given by 

H , (3.2) 

where h^t is a positive constant and tst is the time when the curvature singularity in the star 
appears. When t — )• tgt, H is finite but H and therefore the scalar curvature R = 6H + 12H'^ 
diverge. This implies that in the model in Eq. fl2.13p . the naked curvature singularity is 
generated in the finite future. In the above analysis, we have assumed that the region we 
are considering is almost homogeneous and isotropic. If the assumptions are valid near 
the singularity, the singularity appears simultaneously anywhere in the region. Even if the 
assumptions of the homogeneity and isotropy are broken, the naked singularities appear 
densely in the region. Furthermore, we remark that in general, the matter density becomes 
larger in the region deeper from the surface of the star. Hence, we expect that first, the 
naked curvature singularity could be generated near the center of the star. If the singularity 
generates the attractive force, the star shrinks more, but if the generated force is repulsive, 
there might occur the explosion. If the explosion could occur, however, H must change its 
sign from negative to positive, which seems to be difficult to be realized. 

On the other hand, the model in Eq. fl2.13p could generate the cosmological singularity, 
where the Hubble rate behaves similar to Eq. (13. 2 p as 

H ~ . (3.3) 

In the expanding cosmology, H should be positive and hence the constant hco is positive. 
Moreover, the singularity in the universe occurs at t = tco- The values of tco (and hco) 
in Eq. (13. 3 p and tst (and /igt) in Eq- (13. 2p could be determined by the initial condition for 
H and H. Since we are considering the collapsing star, the absolute values of H and H 
inside the star are expected to be large compared with the values in the expanding universe. 



11 



Therefore, we expect tco > ^st, that is, the curvature singularity in the star occurs before the 
cosmological singularity. 

We note that in order to remove the finite-time future singularity, adding a E? term to 
the form of f{R) works [the first reference in Ref. [l8|]. Similarly, the term can cure 
the curvature singularity under consideration in this paper, as analyzed in Sec. II B. The 
R^ term shifts the value of the potential near the curvature singularity in the scalar tensor 
description. This fact seems to show a close relation between the curvature singularity in 
the star under discussion and the finite-time future singularity in the context of cosmology. 

B. Time scale for the realization of the curvature singularity 

We demonstrate the estimation of the time scale for the realization of the curvature 
singularity by exploring the exponential gravity model in Eq. fl2.24p . In the flat FLRW 
background (13.11) . the gravitational field equations derived from the action in Eq. (12. ip read 



= -\f{R) + ?,(h^ + H) F\R) - 18 (aH^H + HH) F"{R) + /^Vmattcr , (3.4) 
= ^F{R) -[h + 3H^^ F'{R) + 6 (SH^H + 6HH + + F"{R) 

+ 36 UhH + H] ' F"\R) + /t^Pn^atter , (3.5) 



where Pmatter is the pressure of matter. Substituting Eq. (12. 8p with Eq. (I2.24p into Eq. (13.40 . 
we find 

= + /c 1^ + -\-3{h^ + h)^-19, (aH^H + HH^ e't | + «:Vmatter • 

(3.6) 

We examine the case i? — )■ oo. In this case, it follows from Eq. (13. 6 p that H is expressed as 

H^A{t,-t)\n^-^ + H,, (3.7) 

where A is a constant, tg is the time when the curvature singularity appears, ti is a time, 
and Hs is the value of H aX t = tg- When i? — )■ oo, t — )■ and H ^ Hg. In the limit of 
t — J- ts, by combining Eq. (13. 7p with Eq. (13. 6p , we obtain 



= -3H^ + f-[2~ Wt""' ) + -^"-"^ + O (ts - t) , (3.8) 

where we have taken A = Rc/Q and used the fact that the term proportional to H is dominant 
over other terms in Eq. (13.60 because H ^ A/ {t^ — t) and H ~ —A {In [{ts — t) /ti\ + 1}. 
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As an initial condition, we set the initial time t = and H = at t = 0. By deriving the 
solution of = at t = in terms of ti, we have ti = et^. Moreover, by solving Eq. fl3.8p 
in terms of ti, we acquire 



where 



If we regard R 



h = — ^Xe-— +1 , (3.9) 

Itr 



ZK Pmatter 

as the current curvature, we get Rc ~ 12i^^' where = 2.1/ix lO^^^GeV [l 



with h = 0.7 20, 121I] is the current value of the Hubble rate. From ti = etg and Eq. (13. 9p . 
we find 

(3.11) 



where t<, = 2.9 x 10^'' sec 



19| is the age of the flat universe and ^ jt^ = 1.5. 



For the star collapse with the solar mass density Pq = 1.4 gem ^, we find Hg/Hc 
3.9 X 10^'^, where we have used 3H^/k^ = and the critical density pcrit = 3i/^//€^ 



9.2 X 10-28 gcm-3 (l9|. In this case, from Eq. (I3TT|) we obtain 4/4 = 2.9 x lO^^e-i-^^^o^'X. 



If X ~ 0(1), ts/ts, ^ 1, i.e., ts is much smaller than the age of the universe. The reason 
why the time scale for the appearance of the curvature singularity in the exponential model 
is so small originates from the exponential factor e'^^^". We remark that if we apply the 
above analysis to cosmological circumstances, in which Hg ^ He, it follows from Eq. (13. lip 
that ts/tg, ~ 0.27X. As a consequnce, if X ~ 0(1), the time scale for the appearance of the 
curvature singularity in the universe is as long as the age of the universe. 



IV. CONCLUSION 



In the present paper, we have studied a curvature singularity appearing in the star collapse 
for some models of F{R) gravity. We have explored the mechanics to generate the curvature 
singularity. Furthermore, we have explicitly illustrated that the higher derivative term i?" 
(1 < a < 2) could cure the curvature singularity and viable F{R) gravity models could 
become free of such a singularity. It is quite remarkable that same scenario works to cure 
star singularity as well as finite-time future singularity. In addition, we have discussed 
the realization process of the curvature singularity in the viable F{R) gravity model and 
estimated the time scale of its appearance in exponential gravity. It has been shown that in 
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case of the star collapse, the time scale is much shorter than the age of the universe, whereas 
in cosmological circumstances, it is as long as the cosmological time. 
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